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The dynamic spin susceptibility (DSS) has a ubiquitous Lorentzian form in conventional materials 
with weak spin orbit coupling, whose spectral width characterizes the spin relaxation rate. We show 
that DSS has an unusual non-Lorentzian form in topological insulators, which are characterized by 
strong SOC. At zero temperature, the high frequency part of DSS is universal and increases in certain 
directions as cu d_1 with d = 2 and 3 for surface states and Weyl semimetals, respectively, while for 
helical edge states, the interactions renormalize the exponent as d = 2 K — 1 with K the Luttinger- 
liquid parameter. As a result, spin relaxation rate cannot be deduced from the DSS in contrast 
to the case of usual metals, which follows from the strongly entangled spin and charge degrees of 
freedom in these systems. These parallel with the optical conductivity of neutral graphene. 
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Introduction. Strong correlation effects manifest as un¬ 
usual behavior of physical response functions. Of these, 
the frequency and momentum dependent spin suscep¬ 
tibility, x(q,w), played a pivotal role in the study of 
e.g. high-temperature superconductors [TJ, spin-ice com¬ 
pounds Q, and the fundamental description of magnetic 
resonance experiments in correlated systems Q. This re¬ 
sponse function is available experimentally using ac mag¬ 
netization measurements, neutron scattering, magnetic 
resonance, Mossbauer spectroscopy, spin-resolved STM, 
or microwave cavity perturbation experiments. Common 
to these method is that it is difficult to deduce the full 
ui dependent signal, the analysis of experiment therefore 
relies on the theoretical description of the susceptibility. 

The long wavelength spin susceptibility, x(q 0,w), 
called the ac or dynamic spin susceptibility (DSS), in¬ 
dicates dissipative processes and remains in the focus 
of interest when studying the nature of correlations in 
emergent materials, such as e.g. those manifesting the 
spin-liquid phase !4J. DSS is also important in identi¬ 
fying the transition temperature of spin-glasses |tf] and 
superconductors characterizing superparamagnetism 
of small ferromagnetic nanoparticles[7J, or examining the 
nature of magnetic phase transitions. Another highly rel¬ 
evant reason to study DSS is that it provides a measure 
of spin-relaxation rate, whose knowledge is in turn im¬ 
portant for spintronics applications [§j. DSS is character¬ 
ized in the usual materials (where spin-orbit interaction 
is small) by a Lorentzian, which is peaked at the Zeeman 
energy and whose linewidth provides a direct measure of 
the spin-relaxation rate. 

A common feature of the previous cases is the weak 
spin orbit coupling (SOC). However, SOC is usually the 
dominant energy scale in topological insulators[3, |l0( 
which strongly entangles their magnetic properties with 
their charge response. As a result, unusual, non- 
Lorentzian behavior of the DSS might occur. Here, we 


study DSS in topological insulators in the full temper¬ 
ature, doping, magnetic field and frequency range. We 
do find a non-Lorentzian form of the DSS and most sur¬ 
prisingly a non-zero value of the DSS even in the large 
frequency limit. This implies that the spin-relaxation 
rate cannot be determined from the DSS, much as its 
knowledge is desired for prospective spintronics applica¬ 
tions. This result is understood in analogy to the case 
of optical conductivity of neutral graphene: it does not 
follow the usual Drude-Lorentz form due to two-band 
excitations, therefore it cannot be used to determine the 
momentum relaxation rate (llj| . We argue that the non- 
Lorentzian form of response functions is a new hallmark 
of topological insulators. 

In the following, we study four known realizations of 
topological insulators: i) the spin Hall edge state, ii) its 
interacting counterpart, the helical liquid in ID, iii) 2D 
helical Dirac fermions, iv) and the Weyl semimetal in 3D. 

ID Dirac Hamiltonian: the spin-Hall edge state. We 
consider the spin-filtered edge states of a quantum spin- 
Hall insulator[12l-ll4]. whose effective Hamiltonian is 


Hid = vS z p + A S x 


(1) 


with energy spectrum E± ( p ) = ± sj ( vp ) 2 + A 2 and A is 
the Zeeman term from a static magnetic field. 

The DSS requires the calculation of the spin response 
function, which reads in the time domain as 


Xab{t) = l&(t)(S a (t)S b (0) - S b (0)S a (t)), (2) 

where a,b = x, y , or z and S a (t ) = 
exp(—iHidt)S a exp(iHidt) can be calculated using 
the matrix structure of Hid, similarly to Ref. [l5[. A 
given momentum plays the role of an effective magnetic 
field, which acts on the physical spin. Therefore, the 
knowledge of the Xab{t) correlator yields directly the 
DSS. Using the eigenfunctions of Eq. CD. the time 
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dependent correlation function for a given momentum p 
is calculated, yielding the imaginary part of the DSS at 
half filling and T = 0 after Fourier transformation as 
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where F(x) = [(arctan(vT __ x 2 /^) — 7r/2)0(l — |x|) + 
atanh(\/i 2 — l/x)0(|x| — l)]/-y/|i 2 — 1|. The DSS is thus 
strongly anisotropic, contains an off-diagonal term and 
deviates from the ideal Lorentzian form. Depending on 
the geometry, the DSS diverges or vanishes at the gap 
edge and approaches a finite constant value or vanishes 
with increasing frequency. 

In the A = 0 limit, S z is conserved ([S z ,H ld ] = 0), 
therefore x"*M = x" 2 (w) = 0, while x"zM = Xy V (u) = 
sgn(w)/4u, which is the typical density of states in ID. 

The electric current operator is given by j x = evS z , 
therefore the optical conductivity of the spin-Hall edge 
state measure directly x" 2 (w)/w. Note that the other 
components of DSS are not accessible by optical means. 
Additionally, a finite ac electric current can be induced 
along the edge in the presence of an ac magnetic field 
in the y direction due to the finite value of y,, 2 (o;), as 
a manifestation of the magnetoelectric effect [9, 1101 . In 
particular, Rexi /2 (w > 2A) ~ y/~K /(w — 2A). 

The results for the diagonal susceptibilities can be ex¬ 
tended to finite doping and temperature by multiplying 
the calculated Xaa s by sgn(w)(/(^ - f) - f(p+ f)) (ex¬ 
cept for the case of the helical liquid), where f(E ) = 
l/(exp (E/ksT) + 1) with T the temperature and p the 
chemical potential. At T = 0, a finite chemical potential 
introduces an additional gap of 2\p\, and leaves the rest 
intact. At high temperature, it gives a |w|/4T multiplica¬ 
tive factor to the susceptibilities. 

A finite perpendicular magnetic field A opens up a 
gap in the spectrum, and the resulting state becomes im¬ 
mune with respect to interactions as long as \p\ -C A. In 
the absence of the gap, the density of states is finite for 
arbitrary chemical potential, and the interactions pro¬ 
foundly alter the low energy excitations, as is customary 
in 1D[16|. The results obtained below apply also in the 
case of a finite gap, unless y ~ A ~ ui. 

Helical liquid. The helical edge state of the spin-Hall 
insulator forms a helical liquid, when electron-electron 
interaction is taken into account, resembling to a spinless 
Luttinger liquid (LL)[IM3- The Hamiltonian in Eq. 
m is rewritten in second quantized form as [20] 


H 0 



(x)d x Rf{x) 


L+(x)d x L l (x )) , 


(4) 


which is a peculiar half of a spinful LL, lacking the Hq 
and Lf operators. 

The time reversal invariant electron-electron interac¬ 
tion consists of the chiral (<74) and the forward scattering 
( 32 ) terms, 



(5) 


with n-f(x) = R^(x)R^(x) and nj,( x) = L^{x)L^{x). 
These interactions give rise to Luttinger liquid 
behaviour [UH3 with LL parameter K and renormalized 
velocity vf, and K = 1 and vf = v in the non-interacting 
limit. The bosonized Hamiltonian reads as 


H= V -f 

4-7T 
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with the dual fields 9 and </?, satisfying [ip(x),8(y)\ = 
*fsgn {y-x). 



FIG. 1. (Color online) The two possible spin-flip processes 
in the helical liquid, the blue (up spin) and red (down spin) 
lines denote the bare, spin filtered dispersion. The q = 0 
process, corresponding the vertical magenta line, is absent 
in a normal LL and requires a finite frequency threshold 2/x, 
while the green arrow denotes a gapless, q = 2fcp momentum 
transfer process, which does not contribute to DSS, except for 
y. = 0, when these two processes coincide. 


The DSS of the helical liquid is evaluated similarly to 
the 2 kp charge susceptibility of a spinless LL[l6J. The 
spin flip operator is translated to the bosonic language as 
Rf (x)Li(x) ~ exp(—2 ikpx + 2 iip(x)). In the absence of 
perpendicular magnetic field, we obtain x 22 (uO — 0 and 
X"*M = X1M as 
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where B{x , y) = T(x)T(?/)/r(a; + y) is the Euler integral 
of the first kind with T(x) being the Euler’s integral of the 
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second kind [ 23 ], a is a short distance regulator and v F /ct 
represent a high energy cutoff and is shown in Fig. [2] for 
some representative cases. At T = 0, Eq. 0 exhibits 
the typical power law correlation function of a LL as 
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while in the high temperature limit with T ui, /x, it 
yields 
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In spite of the formal similarity to the 2 k F: finite fre¬ 
quency response of normal LLs, Eq. 0 describes a com¬ 
pletely different physical process, which usually involves 
high energy transfer and is beyond the realm of the LL 
paradigm. While the former is gapless in w and accounts 
for a ’’horizontal” interband process with 2 k F momentum 
transfer, the latter stems from a q = 0 ’’vertical” inter¬ 
band transition and is gapped at T = 0 with the thresh¬ 
old frequency of interband transition 2/x, as shown in Fig. 
[Q Only at p = 0, these two processes become identical. 
By the replacement 2p —> 2p ± vq in Eq. 0, the full 
wavevector dependence of the dynamical susceptibility is 
obtained. This indicates that the chemical potential de¬ 
pendent DSS is equivalent to measure the full wavevector 
dependent susceptibility, accessible by e.g. neutron scat¬ 
tering. Similar results were obtained for the dynamical 
density response function as well [ 2 ^. 

This very broad spin response is reminiscent of that 
in the XXZ Heisenberg model[16j, which describes frozen 
charge degrees of freedom due to the strong on site repul¬ 
sion between electrons. The helical liquid, on the other 
hand, operates in the opposite, weakly interacting itiner¬ 
ant electron limit, but the strong SOC entangles the spin 
excitations with the charge degrees of freedom, resulting 
in a broad signal. 

In particular, a strongly repulsive helical liquid with 
K <Cl produces significantly larger spin responses as 
opposed to its weakly or attractively interacting coun¬ 
terpart: the (2-7T aT/ v f ) 2K ~ 2 factor significantly en¬ 
hances/suppresses the spin susceptibility in the repulsive 
(K < l)/attractive (K > 1) case. For K = 1, our pre¬ 
vious expressions for the non-interacting case are recov¬ 
ered. 


Eq. 0 is to be contrasted to the spin response of 
a spinful LL, which in the presence of SU(2) invari¬ 
ant interactions, reduces to uj5(oj ± B ) with B the Zee- 
man field, in spite of the fractionalization of the original 
fermionic excitations into new type of collective bosonic 
modes. Departures from this highly idealized limit im¬ 
ply the inclusion of various SOC terms into the LL 
Hamiltonian [ 2 ^, [2dj| as a weak perturbation on the band 


structure. Our starting point, on the other hand, is the 
completely opposite situation, when the SOC determines 
and dominates the band structure, therefore the SU(2) 
spin rotational symmetry is severely broken and cannot 
be considered as a weak perturbation. 



hj/2 p 


FIG. 2. (Color online) The dynamical spin susceptibility of 
the helical liquid is shown for T = 0 (solid lines) and T = yu/2 
(dashed lines) for several values of the LL parameter. 


2D Dirac Hamiltonian. By increasing the dimension¬ 
ality, the surface states of 3D topological insulators is 
described by the familiar Dirac equation^, given by 

H 2 d=v{S x p y -S y p x ) +AS Z , (10) 


where A is a mass gap, stemming from a thin ferromag¬ 
netic film covering the surface of TI or by a perpen- 
dic ular magnet ic field. The eigenenergies are E±( p) = 
±\/Jyp) 2 TA 2 . 

The time dependent correlation function is obtained 
similarly to the ID case, and the DSS at T = 0 and half 
filling is 
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and Xy.y(w) = y" x (u;). Note that x" y (w) is responsible 
to the ’’half quantum Hall effect”, i.e. the e 2 /2h Hall 
conductivity in topological insulators[9j. Since the elec¬ 
tric current operator is related to the spin due to the 
strong SOC, the in-plane optical conductivity satisfies 
<t(w) ~ y"x( w )/ CiJ ; an d this also agrees with the inter¬ 
band contribution to the optical conductivity of (gapped) 
monolayer graohene(25j. While X xx ,xi /( w ) is measurable 
by optical means as well, the zz component can only be 
probed by magnetic susceptibility measurements. In the 
A = 0 limit, the relation x- z (u;) = 2\ X x{u) = holds 
where the last expression is the typical density of states 
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of e.g. graphene [26}. The factor 2 follows from the spin 
structure of Eq. m- s z sees two perpendicular spin 
components (x and y), which contribute to the response, 
while an in plane component feels only the other in-plane 
component but not S z . Qualitatively similar susceptibil¬ 
ities were derived in Ref. [27 1. 

The effect of a short range electron-electron interaction 
(e.g. Hubbard model) is practically negligible here, as it 
is termed irrelevant in the renormalization group sense 
and can only renormalize the band parameters in the 
weak coupling limit. 



FIG. 3. (Color online) The dynamical spin susceptibility of 
the 2D topological surface state is shown for T = 0 (solid 
lines) and T = 5A (dashed lines) at half filling. 


Weyl semimetal. Inspired by the exciting physics of 
graphene and topological insulators, nodal semimetals in 
3D are currently under investigation [28l - [30j . The Weyl 
Hamiltonian exhausts all three spin operators as 


- V ( S x p x T SyPy -f- S z p z ) . (12) 


The Zeeman energy simply shifts the position of the zero 
energy state in the momentum space and does not open 
a gap in the above Hamiltonian. The DSS follows from 
Eqs. (1111) . after replacing A with k z and performing the 
k z integral, becoming isotropic and diagonal as 


X» 


247TU 3 ’ 


(13) 


being proportional to the density of states of Weyl 
semimetals. Similarly to the previous cases, the opti¬ 
cal conductivity follows as cr(cu) ~ x"(w)/w ~ u> as in 

Ref. [Si. 

Detection. Experimentally, the DSS is directly mea¬ 
sured by the electron spin resonance (ESR) method, 
whose signal intensity is |32|: 


Iaa M = ^WaaMU, ( 14 ) 

where po is the permeability of the vacuum, V is the 
sample volume. Usually, the conventional ESR method 


together with the nuclear magnetic resonance (NMR) in 
solid state systems has limited importance in 2D and es¬ 
pecially ID due to the small number of available states 
(small density of states compared to 3D), which results in 
weak signals. Nevertheless, by considering an ensemble 
of ID nanowires and crystals, the ESR signal can possibly 
be detected similarly to the NMR spectra^ of related 
materials. Additionally, one can also use the recently pro¬ 
posed source-probe setup to measure the DSSjUj. The 
DSS is accessible in a cold atomic realization of these 
states (see e.g. Ref. [35j), featuring also the tunability of 
the interaction strength by standard techniques ( 3 ^], by 
measuring the spin-sensitive Bragg signal, yielding the 
spin-structure factor. 
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